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In the braneworld scenario, the four dimensional effective Einstein equation 

CN \ has extra terms which arise from the embedding of the 3-brane in the bulk. These 

\ non-local effects, generated by the free gravitational field of the bulk, may provide 

an explanation for the dynamics of the neutral hydrogen clouds at large distances 

f-^. \ from the galactic center, which is usually explained by postulating the existence 

of the dark matter. Starting from the observational fact of the constancy of the 
galactic rotation curves in spiral galaxies, we obtain the exact galactic metric in 
the flat rotation curves region in the brane world scenario. The basic physical 
\ parameters describing the non-local effects (dark radiation and dark pressure) 

are obtained, and their behavior is considered in detail. Due to the presence of 
the bulk effects, the flat rotation curves could extend several hundred kpc. The 

"55 ! limiting radius for which bulk effects are important is estimated and compared 

with the numerical values of the truncation parameter of the dark matter halos, 

.^J \ obtained from weak lensing observations. There is a relatively good agreement 

between the predictions of the model and observations. The deflection of photons 
passing through the region where the rotation curves are flat is also considered 
and the bending angle of light is computed. The bending angle predicted by 
the brane world models is much larger than that predicted by standard general 
relativistic and dark matter models. The angular radii of the Einstein rings are 
obtained in the small angles approximation. The predictions of the brane world 
model for the tangential shear are compared with the observational data obtained 
in the weak lensing of galaxies in the Red-Sequence Cluster Survey. Therefore 
the study of the light deflection by galaxies and the gravitational lensing could 
discriminate between the different dynamical laws proposed to model the motion 
of particles at the galactic level and the standard dark matter models. 



Subject headings: galactic rotation curves: gravitational lensing: brane world 
models 
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1. Introduction 

The rotation curves for galaxies or galaxy clusters should, according to Newton's gravi- 
tation theory, show a Keplerian decrease with distance r of the orbital rotational speed Vt g at 
the rim of the luminous matter, v^ g oc M(r)/r, where M(r) is the dynamical mass. However, 
one observes instead rather flat rotation curves (Rubin et al. 1980; Binney and Tremaine 
1987; Persic et al. 1996; Boriello and Salucci 2001). Observations show that the rotational 
velocities increase near the center of the galaxy and then remain nearly constant at a value 
of v t goo ~ 200 — 300 km/s. This leads to a general mass profile M(r) « rvf goo /G (Binney and 
Tremaine 1987). Consequently, the mass within a distance r from the center of the galaxy 
increases linearly with r, even at large distances, where very little luminous matter can be 
detected. 

This behavior of the galactic rotation curves is explained by postulating the existence of 
some dark (invisible) matter, distributed in a spherical halo around the galaxies. The dark 
matter is assumed to be a cold, pressureless medium. There are many possible candidates for 
dark matter, the most popular ones being the weekly interacting massive particles (WIMP) 
(for a recent review of the particle physics aspects of dark matter see Overduin and Wesson 
(2004)). Their interaction cross section with normal baryonic matter, while extremely small, 
are expected to be non-zero and we may expect to detect them directly. It has also been 
suggested that the dark matter in the Universe might be composed of superheavy particles, 
with mass > 10 10 GeV. But observational results show the dark matter can be composed of 
superheavy particles only if these interact weakly with normal matter or if their mass is above 
10 15 GeV (Albuquerque and Baudis 2003). Scalar fields or other long range coherent fields 
coupled to gravity have also intensively been used to model galactic dark matter (Nucamendi 
et al. 2000; Matos and Guzman 2001; Mielke and Schunk 2002; Cabral-Rosetti et al. 2002; 
Lidsey et al. 2002; Mbelek 2004; Matos et al. 2004; Fuchs and Mielke 2004; Lee and Lee 
2004; Mielke and Peralta 2004; Hernandez et al. 2004). 

However, despite more than 20 years of intense experimental and observational effort, up 
to now no non- gravitational evidence for dark matter has ever been found: no direct evidence 
of it and no annihilation radiation from it. Moreover, accelerator and reactor experiments do 
not support the physics (beyond the standard model) on which the dark matter hypothesis 
is based. 

Therefore, it seems that the possibility that Einstein's (and the Newtonian) gravity 
breaks down at the scale of galaxies cannot be excluded a priori. Several theoretical models, 
based on a modification of Newton's law or of general relativity, have been proposed to 
explain the behavior of the galactic rotation curves. A modified gravitational potential of 
the form = — GM [1 + aexp (— r/r )] / (l + a)r, with a = —0.9 and r ps 30 kpc can 
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explain flat rotational curves for most of the galaxies (Sanders 1984, 1986). 

In an other model, called MOND, and proposed by Milgrom (Milgrom 1983; Bekenstein 
and Milgrom 1984; Milgrom 2002, 2003), the Poisson equation for the gravitational potential 
V 2 = AirGp is replaced by an equation of the form V (x) (|V0| /a )] = 4irGp, where a 
is a fixed constant and fJ>(x) a function satisfying the conditions ji (x) = x for x < < 1 
and n(x) = 1 for a; >> I. The force law, giving the acceleration a of a test particle 
becomes a = for » and a = ^a^a^ for « do, where a at is the usual 
Newtonian acceleration. The rotation curves of the galaxies are predicted to be flat, and 
they can be calculated once the distribution of the baryonic matter is known. A relativistic 
MOND inspired theory was developed by Bekenstein (Bekenstein 2004, 2005). In this theory 
gravitation is mediated by a metric, a scalar field and a 4-vector field, all three dynamical. 

Alternative theoretical models to explain the galactic rotation curves have been proposed 
by Mannheim (1993), Mannheim (1997), Moffat and Sokolov (1996), Moffat (2004) and 
Roberts (2004). 

The detection by Brainerd et al. (1996) of weak, tangential distortion of the images 
of cosmologically distant, faint galaxies due to gravitational lensing by foreground galaxies 
has opened a new possibility of testing the alternative theories of gravitation and the dark 
matter hypothesis. Background galaxies are observed to be tangentially aligned around fore- 
ground galaxies because of the latter population's gravitational lensing effect. The angular 
dependence of the shear signal is consistent with the hypothesis that galaxies are domi- 
nated by approximately isothermal haloes, but can also be explained by assuming some 
phenomenological modifications of the effective (Newtonian) gravitational force at large dis- 
tances (Mortlock and Turner 2001). 

Galaxy-galaxy lensing could provide very powerful constraints for the alternative grav- 
itational theories (for a general review of week lensing and its applications see Hoekstra et 
al. (2002)). The measurements rely on averaging over many background sources and the 
signal is only appreciable at large angular separations from the foreground deflectors. Thus, 
in the absence of dark matter, the foreground galaxies can be regarded as simple lenses, and 
the observational data can be used to constrain the deflection law directly. A weak-lensing 
mass reconstruction of the interacting cluster IE 0657-558 has been presented in Clowe et 
al. (2004). 

Weak lensing has been mainly used to discuss and constrain MOND (Gavazzi 2002; 
Hoekstra et al. 2002; Clowe et al. 2004). Gavazzi (2002) used the cluster of galaxies MS2137- 
23, which presents the most constrained lensing configuration of gravitational images ever 
detected in a distant cluster of galaxies to constrain MOND type models. According to this 
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analysis, the MOND model is not compatible with the observations. The need for much 
more baryons in the MOND model than for the dark matter paradigm implies significant 
dynamical differences between these models, which can be explored at very large radial 
distance. Observational evidence against this model was also reported in Clowe et al. (2004). 
For the interacting cluster IE 0657-558 the observed offsets of the lensing mass peaks from 
the X-ray gas directly demonstrate the presence, and dominance, of dark matter in this 
cluster. However, this proof of dark matter existence holds true even under the assumption 
of MOND. Based on the observed gravitational shear-optical light ratios and the mass peak- 
X-ray gas offsets, the "dark matter" component in a MOND regime, that is, the dynamical 
MOND mass M m , related to the dynamical Newtonian mass M d by the relation M m (r) = 



M d /yl + (ao /a) , would have a total mass that is at least equal to the baryonic mass of 
the galactic system. The observed shear and derived mass of the subclump are significantly 
higher than could be produced by an isothermal sphere with a 212 km/s. 

Weak gravitational lensing is a very promising method for the study of the galactic dark 
matter halos. Results of the studies of weak lensing by galaxies have been published recently 
by Hoekstra et al. (2003, 2004). Using weak lensing the flattening of the galactic dark matter 
halos was observed. These results suggest that the dark matter halos are rounder than the 
light distribution, with a halo ellipticity of the order of ehaio ~ 0.33. The average mass 
profile around galaxies was also studied in the framework of two halo models, the truncated 
isothermal sphere model and the Navarro-Frenk- White model. From the point of view of the 
alternative theories of gravitation, the main implications of these results is that spherical 
halos are excluded with 99% confidence, since most of the alternative theories predict an 
almost isotropic weak lensing signal, which is not observed. 

In a series of recent papers (Harko and Mak 2004; Mak and Harko 2004; Harko and 
Mak 2005), several classes of conformally symmetric solutions of the static gravitational 
field equations in the brane world scenario (Randall and Sundrum 1999a,b), in which our 
Universe is identified to a domain wall in a 5-dimensional anti-de Sitter space-time, have 
been obtained (for a review of brane world models see Maartens (2004)). The static vacuum 
gravitational field equations on the brane depend on the generally unknown Weyl stresses in 
the bulk (a higher dimensional space-time, in which our universe is embedded), which can 
be expressed in terms of two functions, called the dark radiation U and the dark pressure 
P (the projections of the Weyl curvature of the bulk, generating non-local brane stresses) 
(Mak and Harko 2004; Maartens 2004; Harko and Mak 2005). Generally, the vacuum field 
equations on the brane can be reduced to a system of two ordinary differential equations, 
which describe all the geometric properties of the vacuum as functions of the dark pressure 
and dark radiation terms (Harko and Mak 2004). In order to close the system a functional 
relation between U and P is necessary. 
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In Harko and Mak (2004), Mak and Harko (2004) and Harko and Mak (2005) the 
solutions of the gravitational field equations on the brane have been obtained by using some 
methods from Lie group theory. As a group of admissible transformations the one-parameter 
group of conformal motions has been considered. The main advantage of imposing geometric 
self-similarity via a group of conformal motions is that this condition also uniquely fixes 
the mathematical form of the dark radiation and dark pressure terms, respectively, which 
describe the non-local effects induced by the gravitational field of the bulk. Thus there is no 
need to impose an arbitrary relation between the dark radiation and the dark pressure. 

As a possible physical application of the solutions of the spherically symmetric grav- 
itational field equations in the vacuum on the brane the behavior of the angular velocity 
v tg of the test particles in stable circular orbits has been considered (Mak and Harko 2004; 
Harko and Mak 2005). The conformal factor ip, together with two constants of integration, 
uniquely determines the rotational velocity of the particle. In the limit of large radial dis- 
tances and for a particular set of values of the integration constants the angular velocity 
tends to a constant value. This behavior is typical for massive particles (hydrogen clouds) 
outside galaxies (Binney and Tremaine 1987; Persic et al. 1996; Boriello and Salucci 2001), 
and is usually explained by postulating the existence of the dark matter. 

Thus, the rotational galactic curves can be naturally explained in brane world models 
without introducing any additional hypothesis. The galaxy is embedded in a modified, 
spherically symmetric geometry, generated by the non-zero contribution of the Weyl tensor 
from the bulk. The extra-terms, which can be described as the dark radiation term U and the 
dark pressure term P, act as a "dark matter" distribution outside the galaxy. The existence 
of the dark radiation term generates an equivalent mass term Mjj, which is linearly increasing 
with the distance, and proportional to the baryonic mass of the galaxy, Mjj(r) ps M B (r/r ). 
The particles moving in this geometry feel the gravitational effects of U, which can be 
expressed in terms of an equivalent mass. Moreover, the limiting value v tgoo of the angular 
velocity can be obtained as function of the total baryonic mass M B and radius tq of the 
galaxy as v tgoo « (2/v^) y/GM B /r Q + (1/12^3) (GM B /r f /2 (Mak and Harko 2004). 

For a galaxy with baryonic mass of the order 1O 9 M and radius of the order of r 70 
kpc, we have Vt g0 o ~ 287 km/s, which is of the same order of magnitude as the observed 
value of the angular velocity of the galactic rotation curves. In the framework of this model 
all the relevant physical parameters (metric tensor components, dark radiation and dark 
pressure terms) can be obtained as functions of the tangential velocity, and hence they can 
be determined observationally. 

However, imposing a group of conformal motions on the brane implies a major restriction 
on the geometrical structure of the space-time. Therefore it would be important to analyze 
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the behavior of the galactic rotation curves in the brane world model without particular 
assumptions on the geometry of the brane. In the framework of a Newtonian approximation 
of the brane metric the possibility that the galactic rotation curves can be explained by the 
presence of the dark radiation and dark pressure was considered in Pal et al. (2005). 

It is the purpose of this paper to consider the geometric properties of the space-time 
at the galactic level on the brane, and to derive the expressions of the physical quantities 
(dark radiation and dark pressure), which determine the dynamic of the particle in circular 
orbit. Under the assumption of spherical symmetry the basic equations describing the static 
gravitational field on the brane depend on two unknown parameters, the dark radiation and 
the dark pressure. As a starting point in our study we adopt the well-established observation 
of the constancy of the galactic rotation curves far-away from the galactic center. This 
property allows the immediate determination of the exact galactic metric on the brane and, 
consequently, of the mathematical form of the dark radiation and dark pressure terms. The 
behavior of other physical parameters (effective potential, angular momentum and energy) 
is also considered. 

As a direct observational application of the obtained results we consider the bending 
of light and the lensing by galaxies in the flat rotation curves region. It has already been 
pointed out that braneworld black holes could have significantly different gravitational lens- 
ing observational signatures as compared to the Schwarzschild black holes (Whisker 2005). 
A general expression for the bending angle in the flat rotation curves region is derived and 
the deflection of light is studied numerically. The deflection angle depends on the tangential 
velocity of the particles in stable circular orbits, the baryonic mass and the radius of the 
galaxy. An analytic expression for the deflection angle in the first order of approximation 
is also obtained. The size of the radius for which the effects of the extra-dimensions are 
important is also derived as a function of the tangential velocity and of the cosmological 
parameters. The explicit expressions of observationally important parameters, like the tan- 
gential shear are presented. Hence the theoretical predictions of the deflection of light in the 
brane world models can be compared with the observations. 

The predictions of the brane world model are compared with the observational data in 
two cases. First, we compare the observational values of the truncation parameter (the size of 
the dark matter halos), obtained in the framework of the truncated isothermal sphere model 
by Hoekstra et al. (2004) with the size of the region for which bulk effect are important on 
the brane. There is a relatively good agreement between these two quantities, the difference 
between prediction and observation being in the range of 15—20%. Secondly, we compare the 
tangential shear in the present model with the observational values of Hoekstra et al. (2004). 
Consistency between brane world models and observations is obtained if the mass-radius 
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ratio of the galactic baryonic matter is of the order of 10~ 4 5 — 10~ 4 . Therefore the bending 
of light could provide a powerful method to distinguish between models which assume that 
dark matter is a form of unknown matter or is a result of a change in the dynamical laws 
that govern the motion of particles. 

The introduction in the past few years of new observational techniques moved cosmology 
into the era of precision science. From the study of the cosmic microwave background (CMB), 
large scale structure of galaxies (LSS) and distant type la supernovae, a new paradigm of 
cosmology has been established. In this new standard model, the geometry of the Universe 
is flat so that fitotai — lj an d the total density is made up of matter, (comprised of baryons 
(Qb = 0.005 and cold dark matter (Qcdm = 0.25) and dark energy with Q\ = 0.7 (Tegmark 
et al. 2004). In particular, current cosmological data provides a very precise bound on 
the physical dark matter density, O dm /i 2 = 0.115 ± 0.012 (Tegmark et al. 2004). This bound 
provides a very strong input on any particle physics or alternative gravitation model for dark 
matter. On the other hand, dark matter is assumed to play a very important role in galaxy 
and large scale structure formation. Since in the model developed in the present paper the 
basic properties of the galactic rotation curves are explained without resorting to dark matter, 
it seems that this would imply that all the dark matter in the Universe could be related 
to extra-dimensional or purely non-standard gravitational effects, as it has been already 
proposed (Cembranos et al. 2003; Maroto 2004). In these models the branons, hypothetical 
particles attached to the scalar fluctuations of the brane, are decoupled from standard model 
matter, thus playing the role of the dark matter. Therefore, a particle component of the 
cosmological dark matter cannot be excluded a priori even in the framework of brane world 
models. Moreover, non-baryonic particles from the standard model or some of its extensions 
could also contribute substantially to the cosmological dynamics. 

The present paper is organized as follows. The gravitational field equations for a static, 
spherically symmetric vacuum brane are written down in Section II. The definition and main 
properties of the observationally important parameters for the study of the galactic rotation 
curves are presented in Section III. In Section IV we derive the galactic metric and the basic 
physical parameters (dark radiation, dark pressure, effective potential etc.) on the brane. 
The deflection of light in the flat rotation curves region is considered in Section V. In Section 
VI we discuss and conclude our results. 



2. The field equations for a static, spherically symmetric vacuum brane 

We start by considering a 5-dimensional space-time (the bulk), with a single 4-dimensional 
brane, on which gravity is confined. The 4— D brane world (^M, is located at a hy- 
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persurface (B (X a ) = 0) in the 5— D bulk space-time (^M, qab) , of which coordinates are 
described by X A , A = 0, 1, 4. The action of the system is given by (Shiromizu et al. 2000) 



S = S, 



bulk 



brane i 



(i; 



where 



and 



Ibulk 



[ 



I (5) 



R + (5) L m + A £ 



2kl " 



d 5 X, 



''brane 



= f * 



.(5), 



Tf K± + Lbrane V') + ^6 



d x, 



(2) 



(3) 



where fcf 



87rG 5 is the 5-D gravitational constant, ^i? and ®L m are the 5— D scalar 
curvature and matter Lagrangian in the bulk, respectively, x^ 1 , /i — 0, 1, 2, 3 are the induced 
4— D coordinates of the brane, K ± is the trace of the extrinsic curvature on either side of 
the brane and Lj, rane (g a ^,ip) is the 4— D Lagrangian, which is given by a generic functional 
of the brane metric g a p and of the matter fields ip. In the following capital Latin indices run 
in the range 0, ...,4, while Greek indices take the values 0, ...,3. A5 and A are the negative 
vacuum energy densities in the bulk and in the brane, respectively. 

The Einstein field equations in the bulk are given by (Shiromizu et al. 2000; Sasaki et 
al. 2000; Maeda et al. 2003) 



(8) G W = kl^T u , ^Tu = -A^g u + 5(B) [-X b ^ 9lJ + T u ] , 



where 



(5) T/ 



IJ 



-2^ g u ^L r 



(4) 



(5) 



'5^g IJ 

is the energy-momentum tensor of bulk matter fields, while T au is the energy- momentum 
tensor localized on the brane and which is defined by 



Tun = 2 ; h g uuLbrane- 



- av 



Sg 



(6) 



The delta function 5 (B) denotes the localization of brane contribution. In the 5-D 
space-time a brane is a fixed point of the Z 2 symmetry. The basic equations on the brane 
are obtained by projection of the variables onto the brane world. The induced 4— D metric 
is 9u 9u ~~ njnj, where nj is the space-like unit vector field normal to the brane 
hypersurface ^M. In the following we assume that all the matter fields, except gravitation, 
are confined to the brane. This implies that ^L m = 0. 
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Assuming a metric of the form ds 2 = (njrij + gi^dx 1 dx J , with nidx 1 = d\ the unit 
normal to the x =constant hypersurfaces and gjj the induced metric on x =constant hy- 
persurfaces, the effective four-dimensional gravitational equations on the brane (the Gauss 
equation), take the form: 

G>z/ = ~A-g fMU + k^T^y + k^S^ u — E^ v , (7) 
where is the local quadratic energy-momentum correction 

S^u = Y2 TT ^ U ~ 4 T ^" Tra + (3T a/3 T a/3 - T 2 ) , (8) 

and E^y is the non-local effect from the free bulk gravitational field, the transmitted pro- 
jection of the bulk Weyl tensor Ciajb, Ejj = CiAjBn A n B , with the property E L j — > 
E^uSiSj as x -> (Shiromizu et al. 2000; Sasaki et al. 2000; Maeda et al. 2003). We 
have also denoted k\ = 8nG, with G the usual four-dimensional gravitational constant. 

The four-dimensional cosmological constant, A, and the four-dimensional gravitational 
coupling constant, /c 4 , are given by A = k\ (A 5 + /c|A 2 /6) /2 and k\ = k 4 \ b /6, respectively. 
In the limit A^ 1 — > we recover standard general relativity. 

The Einstein equation in the bulk and the Codazzi equation also imply the conservation 
of the energy-momentum tensor of the matter on the brane, D V T^ = 0, where D v denotes 
the brane covariant derivative. Moreover, from the contracted Bianchi identities on the brane 
it follows that the projected Weyl tensor should obey the constraint D V E^ = k^DyS^. 

The symmetry properties of E^ imply that in general we can decompose it irreducibly 
with respect to a chosen 4-velocity field as (Maartens 2004) 



E^y — —k 



U (il^Uy + + P^y + 2Q {p Uy) 



(9) 



where k = k^/k^ h^ u = g^y + u^Uy projects orthogonal to the "dark radiation" term U = 
—k 4 E^yU fJ, u u is a scalar, = k 4 h^E a p a spatial vector and P^ v = —k 4 \h( p a h v ^ — \h^yh al3 ~\ E a ( L 
a spatial, symmetric and trace-free tensor. 

In the case of the vacuum state we have p = p = 0, T^ v = and consequently S^ v = 0. 
Therefore the field equations describing a static brane take the form 

iV = -E^y, (10) 

with the trace R of the Ricci tensor R^ v satisfying the condition R = R^ = 4A. is a 
traceless tensor, E£ = 0. 
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In the vacuum case satisfies the constraint P> V E^ = 0. In an inertial frame at any 
point on the brane we have = Sq and h^ v =diag(0, 1, 1, 1). In a static vacuum Q M = 
and the constraint for takes the form (Germani and Maartens 2001) 

l -D,U + UlA^ + D v P iiV + A^P^ = 0, (11) 

where _D M is the projection (orthogonal to w M ) of the covariant derivative and = ^D^u^ 
is the 4-acceleration. 

In the static spherically symmetric case we may chose A^ = A(r)r^ and P^ u = P(r) (r^r u — , 
where A(r) and P(r) (the "dark pressure") are some scalar functions of the radial distance 
r, and r M is a unit radial vector (Dadhich et al. 2000). 

We choose the static spherically symmetric metric on the brane in the standard form 

rfs 2 = _ e »(r) dt 2 + e X(r) dr 2 + r 2 ^2 + sin 2 ^2) (12) 

Then the gravitational field equations and the effective energy-momentum tensor con- 
servation equation in the vacuum take the form (Harko and Mak 2004; Mak and Harko 2004; 
Harko and Mak 2005) 

x (\ A'\ 1 487rG rr . 



12 u'-\' u'\'\ 32nG 



2U + P r(2U + Py y ' 



where we denoted ' = d/dr. 

Eq. (13) can immediately be integrated to give 



_ A _ 1 d 2GMjj (r) A 2 



= 1 - - - — ^ - ^, (17) 
r r 3 

where Ci is an arbitrary constant of integration, and we denoted 

Mu(r) = TJ - / r 2 U(r)dr. (18) 

The function Mu{r) is the gravitational mass corresponding to the dark radiation term 
(the dark mass). For U — and A = the metric coefficient given by Eq. (17) must tend 
to the standard general relativistic Schwarzschild metric coefficient, which gives C\ = 2GM, 
where M = constant is the mass of the gravitating body. 
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3. Stable circular orbits and frequency shifts in static space-times on the 

brane 

The galactic rotation curves provide the most direct method of analyzing the gravita- 
tional field inside a spiral galaxy. The rotation curves have been determined for a great 
number of spiral galaxies. They are obtained by measuring the frequency shifts z of the 
light emitted from stars and from the 21-cm radiation emission from the neutral gas clouds. 
Usually the astronomers report the resulting z in terms of a velocity field v tg (Binney and 
Tremaine 1987; Persic et al. 1996; Boriello and Salucci 2001). 

The starting point in the analysis of the motion of the stars on the brane is to assume, as 
usual, that stars behave like test particles, which follow geodesies of a static and spherically 
symmetric space-time. Next, we consider two observers Oe and O^, with four- velocities u E 
and u^, respectively. Observer Oe corresponds to the light emitter (i. e., to the stars placed 
at a point Pe of the space-time on the brane), and 0^ represent the detector at point P^ 
located far from the emitter and that can be idealized to correspond to "spatial infinity". 

Without loss of generality, we can assume that the stars move in the galactic plane 
9 = n/2, so that u E = (t, r, 0, fi) , where the dot stands for derivation with respect to the 
affine parameter r along the geodesies. In the timelike case r corresponds to the proper 
time. On the other hand, we suppose that the detector is static (i.e., Ooo's four- velocity is 
tangent to the static Killing field d/dt), and in the chosen coordinate system its four- velocity 
is = (i, 0,0,0)^ (Nucamendi et al 2001). 

The motion of a test particle in the gravitational field on the brane can be described by 
the Lagrangian (Mak and Harko 2004) 

where dQ 2 = d6 2 + sin 2 ^^ 2 . From the Lagrange equations it follows that we have two 
constants of motion, the energy E = e u( - r H and the angular momentum / = r 2 (Lake 2004). 
The condition u^u^ = — 1 gives — 1 = —e u ( r H 2 + e x ^r 2 + r 2 2 and, with the use of the 
constants of motion we obtain 

E 2 = e» +x r 2 + e v (^ + lj . (20) 

This equation shows that the radial motion of the particles on a geodesic is the same as 
that of a particle with position dependent mass and with energy E 2 /2 in ordinary Newtonian 
mechanics moving in the effective potential V e ff (r) = e u ^ (l 2 /r 2 + 1). The conditions for 
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circular orbits dV e ff/dr = and r = lead to (Lake 2004) 

I = o r, E 2 = T . 21 

The rotation curves of spiral galaxies are inferred from the red and blue shifts of the 
emitted radiation by stars moving in circular orbits on both sides of the central region. 
The light signal travels on null geodesies with tangent k^. We may restrict k^ to lie in the 
equatorial plane 6 = tt/2 and evaluate the frequency shift for a light signal emitted from Oe 
in circular orbit and detected by O^. The frequency shift associated to the emission and 
detection of the light signal is given by 

z=l-^, (22) 

where uj = —k^Uj, and the index I refers to emission (J = E) or detection (I = oo) at the 
corresponding space-time point (Nucamendi et al 2001; Lake 2004). Two frequency shifts 
corresponding to maximum and minimum values are associated with light propagation in the 
same and opposite direction of motion of the emitter, respectively. Such shifts are frequency 
shifts of a receding or approaching star, respectively. Using the constancy along the geodesic 
of the product of the Killing field d/dt with a geodesic tangent gives the expressions of the 
two shifts as (Nucamendi et al 2001; Lake 2004) 

Z± ~ v/1^772' [3} 

respectively, where exp [v (r)] represents the value of the metric potential at the radius of 
emission r and exp [u^] represents the corresponding value of exp \y (r)] for r — > oo. It is 
convenient to define two other quantities Zd = (z+ — z_) /2 and z\ = (z+ + -2-) /2, given by 

ZD ( r ) = el"—. M]/2_^22=, (24 ) 



^1 -rz//2' 



e [v x -u(r)]/2 

- w = 1 " TT^T!' (25) 

respectively, which can be easily connected to the observations (Nucamendi et al 2001). 
Za and z D satisfy the relation (z A — l) 2 — z 2 D = exp [2 (u D — u(r))], and thus in princi- 
ple exp \y [r)\ can be obtained directly from the observations. This could provide a direct 
observational test of the galactic geometry, and implicitly, of the brane world and other 
extra-dimensional models. 
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The line element on the brane, given by Eq. (12), can be rewritten in terms of the 
spatial components of the velocity, normalized with the speed of light, measured by an 
inertial observer far from the source, as ds 2 = —dt 2 (1 — v 2 ) (Mak and Harko 2004), where 



dr 
dt 



dQ\ 
~dt) 



(26) 



For a stable circular orbit r — 0, and the tangential velocity of the test particle can be 
expressed as 



2 _ /day 
Vt9 ~^\dt) ■ 

In terms of the conserved quantities the angular velocity is given by 

v tg ~ r 2 £2 ■ 



(27) 



(28) 



With the use of Eqs. (21) we obtain 



rv 



J tg 



(29) 



Thus, the rotational velocity of the test body is determined by the metric coefficient 
exp (u) only. 



4. Galactic metric, dark radiation and dark pressure on the brane 

The tangential velocity v tg of stars moving like test particles around the center of a 
galaxy is not directly measurable, but can be inferred from the redshift z^ observed at 
spatial infinity, for which 1 + z^ = exp [(voo — ^) /2] (1 ± v tg ) / yjl — v 2 g . Because of their 
non-relativistic velocities in galaxies bounded by v tg < (4/3) x 10 -3 , we observe v tg ks z^ 
(as the first part of a geometric series), with the consequence that the lapse exp (v) function 
necessarily tends to unity, i. e., e u e u °° / (l — v 2 g ) — > 1. The observations show that at 
distances large enough from the galactic center vt g ~ constant (Binney and Tremaine 1987; 
Persic et al. 1996; Boriello and Salucci 2001). 

In the following we use this observational constraint to reconstruct the metric of a galaxy 
on the brane. The condition 

v 2 = -rv 1 = constant, (30) 
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immediately allows to find the metric tensor component e v in the flat rotation curves region 
on the brane as 



e = 



(31) 



where Rb is a constant of integration. 



By adding Eqs. (14) and (15) and eliminating the dark radiation term U from the 
resulting equation and Eq. (13) gives the following equation satisfied by the unknown metric 
tensor components on the brane: 



, , „ v a 2i/ 2A' i/X' 2 \ 2 , , 
e I i/ + — + — + — ) - — + 4A = 0. 

2 r r 



2 r 2 / r 2 



(32) 



Substituting v given by Eq. (31) leads to a first order linear differential equation satisfied 
by e~ x and which is given by 



d 
dr 
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v. 



K + 1) + 1 1 



2 1 



4A 



v% + 2 



r 6 + v 2 g + 2r vl + 2' 



J tg 



(33) 



From Eq. (33) we obtain e A as 

1 



+ C 2 r~ 2a - 



2A 



« + 2)« (v 2 g + 2)(a + l) 

where C 2 is an arbitrary constant of integration and we denoted 



a = 



'tg 



K + 1) + 1 



v% + 2 



(34) 



(35) 



Eqs. (31) and (34) give the complete galactic metric on the brane. 

Once the metric tensor components are known, the calculation of the dark radiation 
and dark pressure terms is straightforward. Eq. (13) immediately gives the dark radiation 
U on the brane as 



— C/(r) 



v% + 2 



2 (2 - a) , ' 
(a + 1) * 9 



A 



1 



« + 2) a 



1 



- + (2a-l)C 2 r- iQ - 2 . (36) 



From Eqs. (13) and (14) we obtain the dark pressure on the brane as 



96ttG 
kj\b 



P(r) 



3i/ 



+ 



r 



+ 4A, 



(37) 
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or 



967rG 



P{r) = 



vt„ + 2 



v 2 tg + 1 



2(t 



tg 



a 



5< + 1 



{v 2 g + 2) (a + 1) 



A 



5< + 1 



K + 2 ) 



a 



In 2 



+2Q 



2 (v 2 g + 2)a + 5u 



-2a-2 



v% + 2 



(38) 



The angular momentum and the energy of a star moving in the galactic gravitational 
field on the brane are given by 



I = r 



v tg 



and 



E = 



r 

Rb 



1 " v% 



< 1 



1 - v% 



(39) 



(40) 



respectively. For the effective potential describing the Newtonian motion of a test particle 
we find 



r \ 2 < 1 



(41) 



From Eq. (18) it follows that the dark mass associated to the dark radiation term is 
given by 



2GM V (r) 



1 



3 « + 2) 



2 (2 -a) 
(a+T) 



tg 



Ar 3 



K + 2) a 



-2a-l 



(42) 



At distances relatively close to the galactic center, where the effect of the cosmological 
constant can be neglected, but sufficiently high so that the last term in Eq. (42) is also 
negligible small, the dark mass can be approximated as 

v 2 (v 2 + l) V 2 
2GM V (r) ss t9 \ tg >- r w -^r. (43) 

v tg + 2 2 

The dark mass is linearly increasing with the radial distance from the galactic center. 

The behavior of the metric coefficients and of the dark radiation and pressure in the 
solutions we have obtained depends on two arbitrary constants of integration Rb and C<i- 
Their numerical value can be obtained by assuming the continuity of the metric coefficient 
exp (A) across the vacuum boundary of the galaxy. For simplicity we assume that inside the 
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" normal" (baryonic) luminous matter, with density p B , which form a galaxy the non-local 
effects of the Weyl tensor can be neglected. We define the vacuum boundary r of the galaxy 
(which for simplicity is assumed to have spherical symmetry) by the condition p B (r ) ps 0. 

Therefore at the vacuum boundary of the galaxy the metric coefficients are exp (u) = 
1 - 2GM B /r and exp (-A) = 1 - 2GM B /r , where M B = 4tt JJ" p B (r)r 2 dr is the total 
baryonic mass inside the radius r . The continuity of exp (y) through the surface r = r 
gives 

/ 2GM B \~^ , . 

R b = r (l- —^j , (44) 

while the continuity of exp (—A) fixes the integration constant C 2 as 



r 2a 

Co- ° 



< + 2 



/ 2GM B \ , 9 ,1 2A 9 
\ ro J a a + 1 



(45) 



Thus at the galactic level the metric coefficients and the dark radiation and pressure on 
the brane can be obtained in terms of observable quantities. 



5. Light deflection and lensing by galaxies in brane world models 

One of the ways we could in principle test the galactic metric obtained in the previous 
Section would be by studying the light deflection by galaxies, and in particular by studying 
the deflection of photons passing through the region where the rotation curves are flat. Let us 
consider a photon approaching a galaxy from far distances. We will compute the deflection 
by assuming that the metric is given by Eqs. (31) and (34). 

The bending of light on the brane by the galactic gravitational field results in a deflection 
angle A0 given by 

(^) 6raBe = 2|0(r o )-^oo|-7r, (46) 

where 0oo is the incident direction and r is the coordinate radius of the closest approach to 
the center of the galaxy, and generally (Nucamendi et al 2001) 



(r ) - 0c 



PCX 

Jr 



Mr) 

e 2 



e v{r )-u(r) I | _ ^ 



-1/2 



dr 
r 



(47) 



By taking into account the explicit expressions of the metric tensor components in the 
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flat rotation curves region we obtain 

poo 

M - <Poo = I 

J rn 



dr 



1 _)_ n nr -2a 

(«? 9 +2)a + ° 2 ' («? fl +2)(a+l) 



2\ r 2 



(48) 



By introducing a new variable t] = t/tq and taking into account the matching condition 
given by Eq. (65) we can rewrite Eq. (48) as 



M - 0oc = J 



r/ 



^w 9 )_ ^ 1/2 



d?7 



a+1 



(49) 



In standard general relativity the bending angle by a galaxy is given by (A0) GR 



2 |0 (?"o) — 0oo| — Ti" = 4GM B /r . The total deflection angle is given by {A<j)) total 
( A 0W 



(A0) 



brane 



+ 



In order to estimate the magnitude of the bulk contribution to the bending of light on 
the brane and to compare this contribution to the results from standard dark matter models 
we define the parameter 

5 = \? Z ane , (50) 

(A0) DM 

where (A0) DM is the deflection angle as is usually considered in the standard dark matter 
models. Generally, in the dark halo models, a light ray which goes past the halo without 
entering it propagates entirely in a Schwarzschild metric and the light ray is deflected by 
an angle (A<f>) GR . The bending of the light which passes through the dark matter halo is 
determined by the metric inside the halo and this depends on the assumed properties of the 
dark matter. 

We shall compute the parameter 5 in a semi-realistic model for dark matter, in which 
it is assumed that the galaxy (the baryonic matter) is embedded into an isothermal mass 
distribution (the dark matter), with the density varying as p = a^/2nGr 2 , where a v is the 
line of sight velocity dispersion (Binney and Tremaine 1987). In this model it is assumed 
that the mass distribution of the dark matter is spherically symmetric. In fact, if the rotation 
curve is flat, then the mass distribution must be that of the isothermal sphere, for which we 
also have Vt g = V2<J V . The surface density £ of the isothermal sphere is S (r) = a^/2Gr. 
For this dark matter distribution the bending angle of light is constant and is given by 
(A0) DM = 1-kv\ (Blanford and Narayan 1992). 
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Therefore in this model 



(A0) 



brane 



(51) 



The variation of the parameter 5 as a function of the quantity x — ^GM B /r is repre- 
sented in Fig. I. 

In order to find an approximate analytic expression for the deflection of light due to 
the presence of the Weyl stresses on the brane we perform a Taylor series expansion of the 
function F(rj) = (1/ (v 2 g + 2) a + C 2 r]~ 2a - 2Ar 2 r] 2 / (v 2 g + 2) (a + I))' 1 ' 2 near rj = 1. In 
the first order we obtain 



F( V ) 



+ C 2 r 2a - 



K + 2) a 



2Ar 2 



« + 2) (a + 1) 



1/2 2 C 2 r 2a a + 



2Arg 



(^ 2 9 +2)(a+l) 



(17-I) 



+ C 2 r - 



2a 



1 3/2 



«? 9 +2)(a+l) 



■f.. 



(52) 



Substituting this expression into Eq. (49) gives, after evaluating the integral, the first 
order gravitational deflection angle of the light on the brane as 



(A0) 



brane 



= 2 



7T 



1 +c 2 ^ 



[<+ 2 ) a 



(«? 9 +2)(a+l) 



2(l-<) 



( C 2 a + 



2Arg 



«+2)(a+l) 



rf-i *s 



+ c 2 - 



2Arg 



>? 9 +2)(«+l) 



3/2 



- 7T, 



(53) 



where T(z) is the Euler gamma function, defined as T (z) = J °° t z 1 e t dt. 

Alternatively, some very simple expressions for the parameter 5 can be derived in the 
limits x » v tg an d X << v 1 g i respectively. By assuming that v 2 g << 1 and neglecting all 
terms containing the tangential velocity and the cosmological constant, the light deflection 
angle on the brane becomes 



(A0) 



brane 



V-'W-l) (2 - XT 1 )]" 172 dr, 



IT. 



(54) 
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The integral can be calculated exactly, and we obtain 

V2^ K &) +3^2 [{!,!,! },{§,§} ,4\V2 



(A0) 



brane 



XX 



(1-1) 



1/2 



TT, (55) 



where K{m) is the complete elliptic integral of the first kind given by K{m) = j^ 2 (l — m sin 2 9) 1//2 dd 
and p F q (a; b; is the generalized hypergeometric function defined as 



oo 

p F q (a- 6; z) = (ai)* - (%>)* / - (M fc 



fc=0 



A series expansion of Eq. (55) gives 



(A0) 



X 



brane 



(56) 



(57) 



Therefore for the isothermal dark matter model and in the limit of large \ the parameter 
8 tends to the value xl^ v 1 g i $ ~> x/^ v tg- 

Neglecting x an d the cosmological constant in Eq. (49) gives 



(A0) 



brane 



V 



-1 



[fo 2 -l) (2 + < r 1 )]- 1/2 ^ 



— 7T. 



(5? 



A similar calculation as above gives (A0) ferane ~ v^ g /2. Hence in the limit of small x, for 
the isothermal dark matter model 5 can be approximated as 5 ~ 1/47T. These two limiting 
forms of the parameter 5 are consistent with the numerical results presented in Fig. 1. 

If the effect of the cosmological constant can be neglected, the deflection angle of the 
light in the Weyl stresses dominated region around a galaxy is a function of the tangential 
velocity of the test particles in stable circular orbit only. To obtain the total value of the 
bending angle of the light one must also add the usual general relativistic contribution to 
the bending from the baryonic mass of the galaxy. 

Once the light deflection angle is known, one can study the gravitational lensing on the 
brane in the flat rotation curves region. The lensing geometry is illustrated in Fig. 2. 

The light emitted by the source S is deflected by the lens L (a galaxy in our case) and 
reaches the observer O at an angle 9 to the optic axis OL, instead of (3. The lens L is located 
at a distance Dl to the observer and a distance Dls to the source, respectively, while the 
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observer-source distance is D S - r is the impact factor (distance of closest approach) of the 
photon beam. 

The lens equation is given by (Whisker 2005) 

Drq 

tan (5 = tan 9 - [tan 9 + tan ( A0 - 9)\ . (59) 

By assuming that the angle 9 is small, we have tan6> p=* 9 and the lens equation can be 
written as 

« 9 - ^A0. (60) 
Ds 

In the special case of the perfect alignment of the source, lens and observers, j3 — 0, and 
the azimuthal axial symmetry of the problem yields a ring image, the Einstein ring, with 
angular radius 

d (brane) _ ^ A 0. ( 61 ) 

This equation can be expressed in a more familiar form by taking into account that the 
impact parameter r ps D l 9, which gives 



where 9 E is the angular radius of the Einstein ring in the case of standard general relativity, 
9f R) = y/A(D LS /D s D L )GM B . 

In the case of a galaxy with a heavy isothermal dark matter distribution, the Einstein 
radius of the lens formed in perfect alignment is (Blanford and Narayan 1992) 

9 * M) - i^f) tr (63) 

The ratio 5 of the Einstein's rings angular diameters in the brane world models and in 
the isothermal dark galactic halo model is 

Abrane) , » , x 
A — E — \ V) brane (f>A\ 

Of") 2tt< • lMj 

Therefore the ratio of the angular radii of the Einstein rings in the brane world models 
and in the isothermal dark matter model is given by the same parameter 5 which has been 
already introduced in Eq. (50). The variation of the ratio of the Einstein rings in the two 
models is presented in Fig. 1. 
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6. Discussions and final remarks 

The galactic rotation curves continue to pose a challenge to present day physics as 
one would to have a better understanding of some of the intriguing phenomena associated 
with them, like their universality and the very good correlation between the amount of dark 
matter and the luminous matter in the galaxy. To explain these observations models based on 
particle physics in the framework of Newtonian gravity are the most commonly considered. 

In the present paper we have considered and further developed an alternative view to the 
dark matter problem, namely, that the galactic rotation curves can naturally be explained 
in models in which our Universe is a domain wall (a brane) in a multi-dimensional space- 
time (the bulk). The extra-terms in the gravitational field equations on the brane induce 
a supplementary gravitational interaction, which can account for the observed behavior of 
the galactic rotation curves. By using the simple observational fact of the constancy of 
the galactic rotation curves, the galactic metric and the corresponding Weyl stresses (dark 
radiation and dark pressure) can be completely reconstructed. 

The form of the galactic metric we have obtained in the framework of the brane 
world models differs from what would be naively expected, that is, ds 2 = — (1 + 2$) dt 2 + 
(1 + 2$)~ 1 dr 2 + r 2 [d6 2 + sin 2 6d<f) 2 ), with $ representing the Newtonian potential. This 
form is often implicitly assumed (Pal et al. 2005), and the fact that it is not appropriate for 
the region where the rotation curves are flat can lead to significant errors in the estimation 
of the magnitude of some important physical effects, like, for example, the bending of light 
by the galaxies. 

The observations in spiral galaxies usually determine v tg from the redshift zu, so that 
zd ~ Vtg ~ constant. By assuming that exp (i/^) ps 1, and with the use of Eqs. (24) and 
(31), it follows that the results we have obtained are self-consistent if the condition 

(£f>-<r 1/2 *i, («) 

holds. By using the approximation (r/R b )~ Vt3 ps 1 — v 2 g In (r/R b ), valid for \v 2 g In (r/R b )\ « 
1, it follows that the condition is satisfied if \\n(r/R b )\ « v^ 2 . Since v tg ps 10 -3 — 10~ 4 , 
the approximations we have used are self-consistent as long as — 10 6 << In (r/R b ) « 
10 6 , condition which, for the case of the galaxies, does not impose any practically relevant 
constraint. 

As a second consistency condition we require that the timelike circular geodesies on the 
brane be stable. Let r be a circular orbit and consider a perturbation of it of the form 
r = r + 5, where 5 « r (Lake 2004). Taking expansions of V e ff (r) and exp {y + A) about 
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r = r , it follows from Eq. (20) that 

6 + ]-e u{ro)+x{ra) V^ f (r ) 5 = 0. (66) 

The condition for stability of the simple circular orbits requires V"^ (r ) > (Lake 
2004). This gives 

2 (£) r " 4 [-< (1 - 2t£) rl + Z 2 (3 - 5t£ + 2<)] > 0. (67) 

By neglecting the small terms containing powers of the tangential velocity v tg with 
respect to the unity it follows that the circular orbits are stable if their radius satisfy the 
condition r 2 < 3l 2 /v 2 g . 

In the present model there is a very simple relation between the mass of the inter-galactic 
dark radiation M v and the luminous (baryonic) mass M B of the galaxy. By assuming that 
the tangential velocity of particles in circular orbit is approximately given by v 2 g ~ GMb/tq, 
it follows that M v is related to M B via the following simple scaling relation 

Mu(r) « ^-M B . (68) 

The mass of the dark radiation is proportional to the mass of the galaxy and is linearly 
increasing with the distance to the galactic center. 

If we assume that the flat rotation curves extend indefinitely, the resulting space-time is 
not asymptotically flat, but of de Sitter type. This is due to the presence of the cosmo logical 
constant A on the brane. Observationally, the galactic rotation curves remain flat to the 
farthest distances that can be observed. On the other hand there is a simple way to estimate 
an upper bound for the cutoff of the constancy of the tangential velocities. The idea is to 
consider the point at which the decaying density profile of the dark radiation associated to 
the galaxy becomes smaller than the average energy density of the Universe. Let the value 
of the coordinate radius at the point where the two densities are equal be i?™ ax . Then at 
this point (8nG/kj\ b ) U (R^ x ) = (8ttG/c 2 ) 

Punivi where p U nivC 2 is the mean energy density 

of the universe. Hence we obtain 

Din ax 



< + 1 



" \| ^Puniv ~ ^ \^ ~ Vl] A yjv? g {Vl + 1) + 1 



Vtg- (69) 



The mean density of the universe and the value of the cosmological constant can be 
expressed with the help of the density parameters Q = p un i V /Pcrit and f2 A = Ac 2 /3Hq, 
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respectively, where p cri t = 3Hq/8ttG, with H is the Hubble constant, given by H = lOO/i 
km/sec Mpc, 1/2 < h < 1. Therefore 



A numerical evaluation of i?™ ax requires the knowledge of the ratio of the rotation 
velocities in the flat region and of the basic fundamental cosmological parameters. For v tg 
in the range v tg G (10~ 4 , 10" 3 ) and for Q = 1, Q A = 0.7 we obtain R% ax G (315/T 1 , 3150/1" 1 ) 
kpc. On the other hand the measured flat regions are about i?«2x R pt, where R op t is the 
radius encompassing 83% of the total integrated light of the galaxy (Binney and Tremaine 
1987; Persic et al. 1996; Boriello and Salucci 2001). If we take as a typical value R ss 30 kpc, 
then it follows that R < i?™ ax . However, according to our model, the flat rotation curves 
region should extend far beyond the present measured range. 

An alternative estimation of i?™ ax can be obtained from the observational requirement 
that at the cosmological level the energy density of the dark matter represents a fraction 
Vl m ~ 0.3 of the total energy density of the universe Q — 1. Therefore the dark matter 
contribution inside a radius i?™ ax is given by AnVL m (-R™ ax ) Pcrit/3, which gives 



Therefore, by assuming that the dark radiation contribution to the total energy density 
of the Universe is of the order of Vt m m 0.3 we have #^ ax G (388/i _1 , 3881/i" 1 ) kpc for 

v tg e (io- 4 ,io- 3 ). 

The limiting radius at which the effects of the extra-dimensions extend, far away from 
the baryonic matter distribution, is given in the present model by Eq. (70). In the standard 
dark matter models this radius is called the truncation parameter s, and it describes the 
extent of the dark matter halos. Values of the truncation parameter by weak lensing have 
been obtained for several fiducial galaxies by Hoekstra et al. (2004). In the following we 
compare our results with the observational values of s obtained by fitting the observed 
values with the truncated isothermal sphere model, as discussed in some detail in Hoekstra 
et al. (2004). The truncation parameter s is related to Ru by the relation s = Rujl-K (see 
Eq. (4) in Hoekstra et al. (2004) and Eq. (43) in the present paper). Therefore, generally s 
can be obtained from the relation 



nmax tt- 



-1 




(70) 




(71) 




(72) 
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where a is the velocity dispersion, expressed in km/s. Hence the truncation parameter is 
a simple function of the velocity dispersion and of the cosmological parameters only. For a 
velocity dispersion of a = 146 km/s and with fl m = 0.3, Eq. (72) gives s ~ 245/tT 1 kpc, 
while the truncation size obtained observationally in Hoekstra et al. (2004) is s = 213/i 1 
kpc. For a = 110 km/s we obtain s ~ 184/i~ 1 kpc, while o = 136 km/s gives s ~ 228/i _1 
kpc. All these values are consistent with the observational results reported in Hoekstra et al. 
(2004), the error between prediction and observation being of the order of 20%. We have also 
to mention that the observational values of the truncation parameter depend on the scaling 
relation between the velocity dispersion and the fiducial luminosity of the galaxy. Two cases 
have been considered in Hoekstra et al. (2004), the case in which the luminosity L B does 
not evolve with the redshift z and the case in which Lb scales with z as Lb oc (1 + z). 
Depending on the scaling relation slightly different values of the velocity dispersion and 
truncation parameter are obtained. 

The study of the deflection of light (gravitational lensing) in the flat rotation curves 
region can provide a powerful observational tool for discriminating between standard dark 
matter and brane world models. Due to the fixed form of the galactic metric on the brane, 
the light bending angle is a function of the tangential velocity of particles in stable circular 
orbit and the baryonic mass and radius of the galaxy. The specific form of the bending angle 
is determined by the brane galactic metric, and this form is very different as compared to the 
other dark matter models (long-range self-interacting scalar fields, MOND, non-symmetric 
gravity etc.). When x — ^GM B /r « 1, the gravitational light deflection angle is much 
larger than the value predicted by the standard general relativistic approach. Even when 
we compare our results with standard dark matter models, like the isothermal dark matter 
halo model, we still find significant differences in the lensing effect. Therefore the study of 
the gravitational lensing may provide evidence for the existence of the bulk effects on the 
brane. On the other hand, since in this model there is only baryonic matter, all the physical 
properties at the galactic level are determined by the amount of the luminous matter and 
its distribution. 

Generally, the angular position of the source (3 is related to the angular position of its 
image by the lens equation, (3 = 9+a(0). The standard general relativistic deflection law for a 
point mass in general relativity is a (6) = -6 2 E /6 Mortlock and Turner (2001). In order to test 
alternative theories of gravity a more general point mass deflection law, of the form a{9) = 
~ (°e/°) (°o/d + o ) e ~\ f > 1, was introduced by Mortlock and Turner (2001) (see Hoekstra 
et al. (2002) for a more general parametrization). Here 9 is a parameter which can be related 
to the scale r beyond which the physics becomes non-Newtonian. In the brane world model 



the function a can be represented as a (9) = 



(of"') l» 



(r„/GM B ) (A4.\. From the 
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deflection law one can find the tangential shear of the image, 7 tan (9) = [(da/d9) — a(9)/9] /2. 
For the brane world model we obtain 

f e (GR)\ 2 

Ttan (9) = V E Q2 ' (A0) 6rane . (73) 

By using the definition of the parameter 5 we obtain (A0) 6rane = 5 (A(p) DM = 2TrSvf g , 
where, for simplicity, we have adopted again the isothermal sphere model for "standard " 
dark matter. Moreover, we assume that the tangential velocity is related to the baryonic 
mass by the (approximate) Keplerian value, v 2 g ~ GM B /r . Therefore we obtain for the 
tangential shear the following simple expression 



f€ R) ) 



Tt » («) « 2mS v - > . (74) 

This form of the tangential shear is similar to that one resulting from a modification 
of the deflection law proposed by Hoekstra et al. (2002), and which is given by 7 tan (9) = 
(9 out /9 ) {9 2 E /9 2 ), where 9 out and 9 are some ad hoc parameters. The allowable range for the 
values of 9o, obtained from observations, excludes small values. The brane world model fixes 
the values of these parameters as 9 ou t/9o = 2tt5. Since the ratio 9 ou t/9o is a measure of the 
"mass discrepancy", 5 could also be interpreted as a measure of the extra mass generated by 
the dark radiation. The best fit with the observational data is obtained for 9 out = s, where 
s is the truncation parameter. Therefore 5 can also be expressed as 5 = Rjj/9q. 

The mean shear signal can be measured in galaxy-galaxy lensing observations. We com- 
pare the predictions of the brane world model, given by Eq. (74) with the observational data 
obtained in Hoekstra et al. (2004). By fitting an isothermal sphere model to the tangential 
shear at radii smaller than 2' (corresponding to a maximum radius of ~ 350/1" 1 kpc) the 
mean value for the Einstein radius was found to be (9 E ) = 0. 140 ± 0. 009. For this value 
of the Einstein radius the variation of 7 tan calculated for different values of the parameter 
5 is compared with the observational data presented in Hoekstra et al. (2004), in Fig. 3. 
The best fit with the observational data corresponds to large values of 5, 5 ~ 8 — 12. Since 
the velocity dispersion for the considered sample of galaxies is around <rG (128, 150) km/s, 
corresponding to vt g ~ 150 — 210 km/s, it follows that the ratio of the baryonic mass and 
the radius for these galaxies should be of the order of x £ (1-5 x 10 _5 ,5 x 10~ 5 ). If the 
baryonic mass-radius ratio for the galaxies is significantly smaller than this value, then the 
predictions of this simple brane world model could not fit satisfactorily the existing observa- 
tional data. An estimate of the galactic mass-radius ratio can be obtained by approximating 
the tangential value by its Keplerian value, which gives x ~ 8°" 2 - For a = 150 km/s we 
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obtain \ ~ 2 x 10~ 6 . This simple estimate seems to indicate at first sight that the present 
observational data do not favor the brane world interpretation for dark matter. However, we 
have to mention that the estimation of the ratio of the baryonic mass and the radius from 
the tangential velocity could be affected by significant errors. It would be more realistic 
to assume that x ~ 4/c 2 f^ = 8k 2 <7 2 , where k is a factor describing the uncertainty in the 
determination of the ratio of the "normal" mass and radius of a galaxy from the tangential 
velocity alone. If this factor is of the order of 2, k 2, then the brane world model pre- 
diction could become consistent with the observations. Hence to find a definite answer to 
the possibility of the brane world models to correctly describe the weak lensing data a much 
more careful analysis of data is needed. 

Therefore the study of brane world model galaxy-galaxy lensing and of the dark matter 
halos properties could provide strong constraints on the brane world model and on related 
high energy physics models. 

The measurement of the azimuthally averaged tangential shear around galaxies is robust 
against residual systematics, that is, contributions from a constant or gradient residual shear 
cancel. However, this is not the case for the quadrupole signal (Hoekstra et al. 2004). If 
the lens galaxies are oriented randomly with respect to the residual shear, then the average 
over many lenses will cancel the contribution from systematics. But in real observations 
the uncorrected shapes of the lenses are aligned with the systematic signal. An imperfect 
correction can give rise to a small quadrupole signal, which can mimic a cosmic shear signal. 
Generally, the residual shear has an amplitude 7 and is located with respect to the major 
axis of the lens at an angle 0. The tangential shear 7°^ observed at a point (r, 9) is the 
sum of the lensing signal 7'^ and the contribution from the systematics, 7 tan , so that 
7t°an = liTn + Ttan- 7tan is given by 7tan = -7 cos [2 (0 - 0)] (Hoekstra et al. 2004). One way 
to estimate the flattening of the halo is to measure the shears 7^ at 6 = and ir and 7^ 
at 6 = 7r/2 and 3tt/2, respectively. The observed ratio is f b. s = 7tan/7tan can be written 
as f b s = 7 t ^ +7 cos (20) / 7^ — 7 cos (20) j (Hoekstra et al. 2004). In the framework 
of the simple brane world model we have considered all the shear parameters (tangential as 
well as residual) are proportional to the values corresponding to the standard dark matter 
case. On the other hand the mean value of cos (20) is proportional to the measure a of the 
correlation between the position angle of the lens and the direction of the systematic shear of 
the background galaxies, which is a very small quantity. The brane world model correction to 
a cannot lead to a significant increase to this quantity, due to the large separation distances 
between galaxies. Therefore we expect that f Q b s is an invariant and brane world effects do 
not affect the robustness of the measurement of the average halo shape. 

With the use of the approximate relation (1 + x/n) n ~ exp (x), which is valid for large 
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n, we can write Eq. (44) in the form 

R b ^r exp (xlH g ) • (75) 

In the case of a very massive object with x ~ 1 an d for v tg G (1CT 4 , 10~ 3 ) it follows 
that Rb » r . But for very small values of x, of the same order as v^ g , Rb ~ r . Therefore 
flat rotation curves are specific for particles moving in stable circular orbits around galaxies 
(having x << 1)) while the same phenomenon cannot be detected for very compact (stellar 
or black hole type) objects. 

The measurement of the anisotropy in the lensing signal around galaxies and the de- 
tection of the flattening of the dark matter halos could pose some serious challenges to the 
alternative theories of gravitation (Hoekstra et al. 2002, 2004). In many theories proposing 
modifications of the dynamic law for gravitation the lensing signal caused by the intrinsic 
shapes of the galaxies decreases as oc r -2 and hence it is negligible at large distances. There- 
fore dark halos cannot be modelled as spherical systems. This important observational result 
can clearly discriminate between the gravitational explanations proposed as substitutes for 
dark matter. However, in the framework of the model discussed in the present paper, the 
flattening of the dark matter halos and the anisotropic signal can be easily explained, at 
least at a simple qualitative level. We have obtained all our results under the (unrealistic) 
assumption of the spherical symmetry, with the galaxy (consisting of baryonic matter only) 
modelled as a self-gravitating relativistic sphere. A much more realistic model for a galaxy 
is its representation as a self-gravitating rotating stationary axisymmetric disc, in which 
generally the metric coefficients depend on both polar coordinates r and 9. The brane world 
metric outside this baryonic matter distribution is also anisotropic, and we would naturally 
expect a flattening of the region in which the bulk effects are important, associated with an 
anisotropic (angular dependent) weak lensing signal. The "dark matter" is, in this model, a 
projection of the geometry of the higher dimensional space time far away from the baryonic 
matter distribution, which has to match the geometry of the galaxy. Therefore disc type or 
spheroidal distributions will automatically create flattened and anisotropic (non-spherical) 
geometrical effects. For example, already in standard general relativity in the weak field limit 
a cold disc has the associated gravitational potential — v = $ = — v^ g \n[r (1 + |cos#|) /D], 
where D is a fiducial length scale. The bending of light is anisotropic (polar angle depen- 
dent), as is the focussing effect of the disc (Cai and Shu 2002). Therefore, more realistic 
geometric distributions of the baryonic matter inside the galaxy and of the exterior space- 
time geometry of the brane can easily explain the flattening effect and the anisotropic lensing 
observed for dark matter halos. 

One of the main advantages of the brane world models, as compared to the other alter- 
native explanations of the galactic rotation curves and of the dark matter, is that it can give a 
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systematic and coherent description of the Universe from galactic to cosmological scales. On 
the other hand, in the present model all the relevant physical quantities, including the dark 
energy and the dark pressure, which describe the non-local effects due to the gravitational 
field of the bulk, are expressed in terms of observable parameters (the baryonic mass, the 
radius of the galaxy and the observed flat rotational velocity curves). Therefore this opens 
the possibility of testing the brane world models by using astronomical and astrophysical 
observations at the galactic scale. In this paper we have provided some basic theoretical 
tools necessary for the in depth comparison of the predictions of the brane world model and 
of the observational results. 
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log(X) 

Fig. 1. — The ratio 5 = (A0) femne /2irVt g of the light deflection angle in the galactic metric 
on the brane and the deflection angle in the semi-realistic isothermal galactic dark matter 
model as a function of the parameter x — 4:GM B /r for different values of the tangential 
velocity: v tg = 150 km/s (solid curve), v tg = 210 km/s (dotted curve), v tg = 300 km/s (short 
dashed curve) and v tg = 450 km/s (long dashed curve). For the impact parameter r and 
the cosmological constant A we have adopted the values r = 30 kpc and A = 3 x 10~ 56 
cm -2 , respectively. 
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Fig. 2. — The lensing geometry, showing the location of the observer O, the lensing galaxy 
L and the source S. The deflection angle is A<p. The angular diameter distances D L , D L $ 
and D s are also indicated. 



-33- 




40 60 80 
r [arcseconds] 



Fig. 3. — Comparison of the variation of the tangential shear as a function of the distance 
in the brane world model and the observational data of Hoekstra et al. (2004) (represented 
by points) for different values of the parameter 5: 5 = 6 (solid curve), 5 = 8 (dotted curve), 
5 = 10 (dashed curve), 5 = 12 (long dashed curve). The Einstein radius is 9 E = 0."l4. 



